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ABSTRACT 

A new method for modelling spherically symmetric inhomogeneities is applied to the 
formation of clusters in an expanding Universe. We impose simple initial velocity and 
density perturbations of finite extent and we investigate the subsequent evolution 
of the density field. Photon paths are also calculated, allowing a detailed considera- 
tion of gravitational lensing effects and microwave background anisotropies induced 
by the cluster. We apply the method to modelling high-redshift clusters and, in par- 
ticular, we consider the reported microwave decrement observed towards the quasar 
pair PC1643-I-4631 A&B. We also consider the effect on the primordial microwave 
background power spectrum due to gravitational lensing by a population of massive 
high-redshift clusters. 

Key words: Gravitation - cosmology: theory - cosmology: gravitational lensing - 
cosmic microwave background - quasars: individual: PC1643-I-4631 A&B - galaxies: 
clustering 



1 INTRODUCTION 

■ Several methods have been proposed for modelling the for- 
' mation of galaxy clusters in an expanding Universe. Ide- 
, ally, one would like to model the evolution of arbitrary 
' matter perturbations in a fully general-relativistic manner. 

Unfortunately, this has not proved possible owing to the 
. complication of the calculations involved, and models of 
' the formation of clusters of arbitrary shape employ an ap- 
proximate linearised approach to incorporate the effects of 
gravity (Martmez-Gonzalez, Sanz & Silk 1990; Martmez- 
Gonzalez & Sanz 1990; Chodorowski 1991). Nevertheless, 
many galaxy clusters can be reasonably approximated as el- 
lipsoids, and some examples, such as the Coma cluster, are 
quasi-spherical. This naturally leads to the use of spherical 
symmetry as an approximation, which greatly reduces the 
computational complexity and allows the full incorporation 
of general-relativistic effects. Early approaches to this prob- 
lem were based on the 'Swiss Cheese' model (Rees & Sciama 
1968; Dyer 1976; Kaiser 1982; Nottale 1982 & 1984), whereas 
more recent attempts have used the continuous Tolman- 
Bondi solution (Panek 1992; Arnau et al. 1993; Saez, Arnau 
& FuUana 1993; Arnau, Fullana & Saez 1994; Fullana, Saez 
& Arnau 1994; Quilis & Saez 1998). 

The Swiss Cheese and Tolman-Bondi approaches both 
assume that the cosmological fluid is pressureless. This as- 
sumption may be reasonable for very large structures in 
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which the baryon content is negligible compared to the dark 
matter. For typical clusters, however, the assumption is per- 
haps questionable since 10-30 percent of the mass of the 
cluster may comprise a hot baryon component, with a non- 
negligible pressure, that is gravitationally coupled to the 
dark matter. Nevertheless, calculations by Quilis, Ibafiez & 
Saez (1995), that include a hot gas component, show that 
the eff'ects of pressure are in fact negligible, and that the 
collapse of the cluster is well approximated by assuming a 
pressureless fluid. However, the combined use of spherical 
symmetry and pressureless assumptions lead inevitably to 
the unrealistic situation where the cluster collapses to form 
a singularity. Therefore, for the applications described in this 
paper, we demand the density profile to be realistic at the 
time photons observed today are traversing the cluster. As 
pointed out by Quilis et al. (1995), the assumption of spher- 
ical symmetry introduces another difficulty. As a result of 
the radial structure of the velocity field, gravitational forces 
generate infalling motions that are too rapid and lead to 
fast evolution of the cluster, even when a hot gas compo- 
nent is included. An important consequence of this effect 
is that estimates of cosmic microwave background (CMB) 
anisotropies due to the cluster are overestimated and should 
be regarded as upper limits. However, for high redshift clus- 
ters still in the process of formation, such infalling motions 
may be less unrealistic (Navarro, Frenk & White 1996). 

With these points in mind, in this paper we apply a 
new method for modelling the evolution of spherically sym- 
metric anisotropies in the Universe, assuming a pressure- 
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less fluid. This method is discussed by Lasenby et al. (1998) 
(hereinafter Paper I) and is based on a new, gauge-theoretic 
approach to gravity (Lasenby, Doran & Gull 1998) that pro- 
vides a simple framework in which to investigate spherical 
cluster formation. A brief outline of the method is given in 
Section ^ below. Our intention is to focus on cluster for- 
mation at high redshifts. However, in order to compare our 
results with existing work we first consider low redshift clus- 
ters in some detail. In Sections ^ and ^, we model the for- 
mation of low-redshift clusters with typical characteristics 
and calculate the CMB anisotropies expected from these 
structures. Our results are compared with previous stud- 
ies. In Section we consider the formation of high-redshift 
clusters and discuss their effects on the CMB and the grav- 
itational lensing of background sources. As an application, 
we model a cluster with a view to explaining the microwave 
decrement reported towards the quasar pair PC1643-I-4631 
A&B (Jones et al. 1997). Finally, in Section ^ we consider 
the effects on the power spectrum of primordial CMB fiuc- 
tuations of gravitational lensing by a population of massive 
high-redshift clusters. 



2 THEORETICAL MODEL 

The theoretical model used in this paper is discussed in de- 
tail in Paper I and provides an exact general-relativistic solu- 
tion describing the evolution of a spherically-symmetric dis- 
tribution of pressureless fluid. The 'Newtonian' gauge used 
for describing this physical system is global, having a single 
time coordinate t and a radial coordinate r. The time co- 
ordinate t measures the proper time for observers comoving 
with the fluid. In inhomogeneous regions the radial coordi- 
nate r is related to the strength of the tidal force defined 
by the Riemann tensor. Throughout this paper we employ 
natural units G — c = h — 1, unless stated otherwise. 

The main dynamical variables in this model are the fluid 
density p{t,r), its radial velocity field u{t,r) and a gener- 
alised 'boost' factor r{t,r) given by 



2M{t,r)/r + u{t, 



(1) 



which in the case M{t, r) = reduces to the square of the 
standard special relativistic 7-factor. As explained in Pa- 
per I, M{t, r) is the the total gravitational mass within a 
coordinate radius r and is defined by 



(2) 



M{t,r) = / Atvs p{t,s)ds. 



Given the density and radial velocity fields p{ti,r) and 
u{ti,r) defined at some initial time ti, we may calculate 
r{ti, r) and M{ti,r), which are then conserved along stream- 
lines. The evolution of the fluid density and velocity are then 
easily calculated. 

We choose to impose a very simple finite initial pertur- 
bation in the radial velocity field u{ti,r) as discussed below. 
As explained in Paper I, assuming that the initial velocity 
and density profile arose from primordial perturbations, we 
infer the corresponding initial density profile using the rela- 
tion 



u{ti, r) 



2r 



2H{U 



M{U,r) 



(3) 




Figure 1. Geometrical arrangement defining the variables used 
in the theoretical model. O denotes the observer, S the source 
emitting the photon and C the cluster, which is located at the 
origin of the coordinate system. Both the observer and the source 
are comoving radially with the cosmological fluid. The vector p 
represents the covariant photon momentum. 



which, by differentiating (^), gives 



8tt 



r dr 



(4) 



We note that (^ and (^) are valid for 0,{ti) ~ 1, which will 
always be true if ti is soon after infiation. An important 
respect in which our approach differs from that of Panek 
(1992) or Quilis et al. (1995) is that the density enhance- 
ment resulting from the initial velocity perturbation is au- 
tomatically compensated within a finite region by a slightly 
under-dense region surrounding it. Thus the external Uni- 
verse never feels a gravitational influence from the perturbed 
region. In addition, as explained in Paper I, imposing an 
initial velocity perturbation and working in a fully com- 
pensated manner avoid problems of streamline crossing (i.e. 
shocks) . 

The initial velocity field u{ti,r) is controlled by four 
parameters Ri, a, m and H{ti) as follows. The initial lin- 
ear extent of the perturbed region is Ri and the magnitude 
of the perturbation is controlled by the velocity gradient at 
the origin, which we denote by a. The velocity perturbation 
u{ti, r) for Q <r < Ri is defined by a (2m -I- l)-degree poly- 
nomial of which the first m derivatives are matched at the 
boundaries (i.e. at r = and r — Ri). The fourth parameter 
H{ti) is related to the external Universe. Outside the per- 
turbed region (r > Ri), the fiuid is described by a uniformly 
expanding FRW model. This is straightforwardly achieved 
by setting u{ti,r) — H{ti)r, where H{ti) is the Hubble pa- 
rameter at the initial time ti. As shown in Paper I, choosing 
the velocity gradient at the origin smaller than H{ti) results 
in the region r < Ri moving inwards relative to the Hubble 
fiow and eventually collapsing to form a cluster centred at 
the origin. Conversely, choosing a gradient larger than H{ti) 
would lead to the formation of a void, and the investigation 
of the evolution of such structures is also available in this 
approach. 

Once Ri, a, m and H{ti) are fixed, no further parame- 
ter is required to constrain the initial velocity perturbation 
and the density perturbation is then fully defined by equa- 
tion (^). The initial central density is therefore given by 

3// (to 



piU,0) = 



87r 



■ {4H{t,) - 3a) ^ 



(5) 



and p{ti,r) is represented by a (2m)-degree polynomial for 
< r < Ri. We note that, at r = 0, the fluid evolves 
as a closed FRW Universe with initial density p{ti,0) and 
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velocity gradient a. Outside the perturbed region, for r > 
Ri, the density is uniform and equal to 



Pi = 



3H\ti 
Stt 



(6) 



We note that m has to be greater than 1 in order to obtain 
a sensible compensated density perturbation with a (2m)- 
degree polynomial. 

In addition to the evolution of the fluid, the calculation 
of photon trajectories and redshifts is also straightforward. 
We can parameterise the photon path in terms of the time 
parameter t, so that the path is deflned by r(t) and (j){t 
The geometrical conflguration adopted is illustrated in Fig. 
where x is the angle between the photon's position vector 
and its covariant momentum p. As explained in Paper I, x is 
an observable quantity and is also equal to the angle between 
the photon trajectory and the cluster centre, as measured by 
observers comoving with the fluid. Given an initial position 
{rQ,4>o) and direction xo, the set of first-order differential 
equations given in Paper I may be integrated numerically to 
obtain the subsequent photon path and frequency. 




0.05 0.1 0.15 

R, (h(,"' Mpc) 

Figure 2. The cluster core radius Rc at z = 0.09 as a function 
of the initial perturbation size Ri for various values of the integer 
parameter m (the latter controls the polynomial order - see text). 
The points of intersection between the dashed line and the solid 
lines indicate the suitable value of Ri for each case. 



3 CLUSTER FORMATION 
3.1 Initial conditions 

Using the model outlined above, it is straightforward to in- 
vestigate the formation of spherical clusters. We begin spec- 
ifying the initial conditions such that the resulting cluster 
is similar to those observed. First, we must choose the time 
ti at which the initial perturbation is applied. Throughout 
this paper we take ti to represent the epoch z = 10^ and 
we assume that the Universe is at critical density so that 
ri(t) — 1 for all values of t. Unless stated otherwise, we per- 
form simulations for H{to) — 100 km s~^ Mpc~^ where to is 
the present time. Results will be quoted using the reduced 
Hubble parameter h{t) which is equal to H{t) expressed in 
units of 100 km s~^ Mpc~^. Furthermore, quantities evalu- 
ated at t = to will be written simply with a zero subscript, 
e.g. f2(to) = f2o, h{to) = ho, and so on. The free parameters 
of our model are H{ti), which defines the external Universe, 
and a, Ri and m, which d eterm i ne th e initial perturbation. 
The following paragraphs (B.1.1, 3.1.2 & 3.1.3) describe how 
we fix these parameters. 



3.1.1 External Universe 

1 and ho, the parameter H{ti) is 



Given Q,o ~ 
flned. For t 
lO'^/io km s 



representing z 
^ Mpc-^ 



10^ we find H{t,) 



fully de- 
= 3.167 X 



3.1.2 Final appearance of the cluster 

The only guide to fixing the initial parameters a, Ri and 
m is the final properties of the cluster, as observed today. 
In this section, we model a very rich Abell cluster, in order 
to compare our results with those obtained by Quilis et al. 
(1995) and Panek (1992). As our standard configuration, we 
consider a cluster at a redshift z — 0.09, with a core ra- 
dius Rc = 0.23/io ^ Mpc (where Rc is defined as the radius 
at which the cluster density falls to one-half its maximum 




Figure 3. The total gravitational cluster mass at 2 = 0.09 within 
a sphere of radius r in unit of 10^^ h^^ Mq , as a function of r for 
various values of m. Ri has been fixed using results from Figure El 



value). A redshift z = 0.09 corresponds to an effective dis- 
tance of D = 250/Iq ^ Mpc. The central baryonic number 



density of the cluster is taken to be Tie 



.,„4t1/2 -3 

10 ftn m . 



We 



further assume that, at all times, the baryon component con- 
tributes lOhg^^^ per cent of the total mass, the remainder 
being dark matter. Finally, because we are interested in rich 
clusters, the total gravitational mass Mc within r < J?c is 
taken to be greater than 10^^ h^^ Mq. Note that the values 
of Rc, Uc and Ale are defined at the epoch when a photon 
observed today was traversing the cluster. 

We should discuss our choices regarding the h- 
dependences stated above. The /i-dependence of the core 
radius ensures that the observed angular size of the clus- 
ter is independent of the Hubble parameter and equal here 
to 3.4 arcmin. The /i-dependence of the central baryonic 
number density is chosen so that the cluster's X-ray fiux 
is also independent of ho (Peebles 1993; Jones & Forman 
1984). Hence the three observational properties of the clus- 
ter, namely its redshift, angular size and X-ray flux do not 
depend on our choice of ho. The /i-dependence of the dark 
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Figure 4. Initial relative velocity u(ti,r)—H{ti)r for the m = 2, 3 
and 4 models at ^ = 10''. 




m=4 m=3 
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Figure 5. Initial density profiles p{ti,r) for the m = 2, 3 and 4 



models at 



10'=' 



matter ratio gives a total gravitational mass which scales as 
as expected within the isothermal sphere approximation 
(e.g. Peebles 1993). In this case the tot al d ensity follows the 



power law form p oc 1/r (see Section 3.3) and the baryon 



mass alone scales as h. 



-5/2 



We note that the distance D = 250/1(7^ Mpc is some- 
what larger than the value of lOOhg^ Mpc used by Quilis et 
al. (1995). This is because our simulations show that in the 
case of the cluster described above, but located at a distance 
of 100/Iq ^ Mpc, the observer is marginally decoupled from 
the Hubble flow (particularly for the m — 2 model). Thus 
in order that the observer resides in the external Universe it 
was necessary to place the cluster at a larger distance (see 
Section 13). 



3.1.3 Fixing the initial perturbation 

At the origin (r = 0), the fluid evolves like a closed FRW 
Universe with initial density p{ti,Q) given in ^ and ve- 
locity gradient a. Therefore a is directly constrained by re- 
al z = 0.09. We find that 



m 



Rc = 0.23/io ^ Mpc and Mc > 10' 



Mp) at z = 0.09, 



quiring that ric = l(fh^^'^ 
1 — [a/H{ti)\ ~ 6 X 10~^. The further requirements, namely 

that allow us to fix both Ri and m as follows. We find 

0.23/1(7^ Mpc, Ri 
0.057/io\ 0.044/io^ 



(see Fig. ^ that in order to satisfy i?c 
should be equal to G.21h^^ , 0.089^0" 

and 0.036/1(7^ Mpc for m = 2, 3, 4, 5 and 6, respectively. 
Because we require that Mc > IQ^^h^^ M0 at z — 0.09, we 
can see from Fig. ^that only m = 2, m = 3 and m = 4 cases 
are suitable. The initial velocity and density perturbations 
for the m = 2, 3 and 4 models are shown in Figs. ^ and ^ re- 
spectively. Each of these perturbation give rise to the same 
cluster properties at z = 0.09, namely: Uc = 10^/iJ''^ rn"'', 
i?c = 0.23h~^ Mpc and Mc > lO^^h'^ Mq. We note that 
the m = 2 model requires a much larger perturbation. This 
could be explained by the fact that the initial density profile 
for larger m is flatter at the origin, concentrating more mass 
at the centre and therefore allowing for a faster collapse. As 
a typical model, most of the figures presented in this p aper 



are for m = 3 while all the numerical results of Sections 3.2 
3.3 and W will be quoted for the three cases m = 2, 3 and 4 



3.2 Fluid evolution 

In order to appreciate the results presented in this paper, 
in particular the effect on the CMB induced by the evolving 
cluster, we present a series of figures (Figs. ^ ^ ^ and ^ 
illustrating various fluid quantities as measured by a set of 
observers comoving with the fluid, on the path of a photon 
traversing the centre of the cluster. These figures allow a 
check that our model is behaving sensibly. 

Fig. 1^ shows the baryonic number density, as a function 
of cosmic time t, along the path of a photon which travels 
straight through the centre of the cluster and reaches the 
observer at the current epoch to; the curve displayed is for 
m = 3. For convenience we take the current epoch to — 0, 
so that all times are expressed in Myr prior to today. From 
the figure we see that the photon traversed the cluster in 
approximately 5 Myr, encountering a maximum baryonic 
number density of 10'* ftp ^''^ m"'' at a time tc — 738/i(7^ Myr 
ago. 

Similarly, Figs. ^ & ^ show respectively the fluid ra- 
dial velocity field u{t,r) and velocity gradient du/dr along 
the photon path for m — 3. Several features in these plots 
should be noted. It is clear from Fig. ^ that from about 
500/i(7* Myr prior to the current epoch, the velocity gradi- 
ent is merely that due to the slowing expansion of the ex- 
ternal Universe, and so the photon is free-streaming towards 
the observer. As required, this velocity gradient tends to the 
value Ho = lOO/io km s"^ Mpc"* at t = 0. From Fig. 0, we 
see that, even though the photon experiences a perturbation 
due to the cluster up to 500hg * Myr prior to today, it is only 
when in the inner part of the cluster, between SlOh^^ and 
675/i(7* Myr ago, that the fluid velocity is directed inwards. 
In the outer parts of the perturbation the fluid is still moving 
outwards, although it is, of course, collapsing with respect 
to the Hubble flow. The fluid velocity relative to the Hubble 
flow is shown in Fig. ^ 

For models with m 7^ 3, the corresponding results are 
very similar to those shown in Figs. ^, ^ ^ and ^ but with 
a larger or smaller perturbed region for the models m — 2 
and m = 4 respectively. 
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Figure 6. The baryonic number density, as a function of cosmic 
time t, along a photon path traversing the centre of the cluster. 
The time t = corresponds to the current epoch and we assume 
m = 3. 
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Figure 7. The radial velocity field u{t, r) along the photon tra- 
jectory. The time i = corresponds to the current epoch and we 
assume m = 3. 
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Figure 8. As in Fig. but for the velocity gradient du/dr along 
the photon path. 
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Figure 9. As in Fig. ^ but for the relative velocity u(t, r) — H{t)r 
along the photon path. 



3.3 Cluster characteristics 

Before giving numerical values of tlie cluster cliaracteristics 
and comparing them witii previous worlc, we first clieclc tliat 
the obtained cluster profile has a realistic shape. We consider 
here a broad family of density profiles (Zhao 1996, Kravtsov 
et al. 1998) given by 



n(r) — 



no 



(7) 



We note that this family includes the three most widely 
used types of density distribution: (i) Profiles derived from 
A'^-body codes and proposed by Navarro et al. (1996), which 
diverge at the origin, hereafter NFW profiles. These corre- 
spond to (a,/3, 7) = (1, 1, 1) in equation (^). For this case, 
no/4 is the density at r = _R; (ii) Isothermal profiles, also 
known as King profiles (Rood et al. 1972, Sarazin 1988), 
which have a finite central density n = no and for which 
(a,/3,7) — (2,2/3,0). Here, the density profile follows the 
very simple form 



n(r) — 



no 



1 



(8) 



and R is the radius at which n = no/2; (iii) We define here 
/3-models as being a generalisation of King profiles. They are 
widely used to fit observed profiles and rotation curves and 
correspond to (a,/3, 7) — (2,/3, 0). 

In order to compare our simulated cluster profiles with 
equation (Q), we consider the variation with radius r of 
the number density nc(r) of our clusters at the time tc = 
738/Iq"^ Myr ago, when the photon reaches the point of 
maximum density along its trajectory. This profile is plot- 
ted as the solid line in Fig. hd for the case m = 3. Per- 
forming a simple least-squares fit for the three parameters 
of the /3-models from r = to r = l.bhg^ Mpc, we find 
no = 1.04 ± 0.08 X h^^^'^lO'^ m"^ R = 0.25 ± OMh^^ Mpc 
and P — 0.75 ± 0.15, where the quoted errors are one-sigma 
limits. The best-fit /3-model is shown as the dashed line in 
Fig. |l^. The m — 4 model gives rise to a density profile 
somewhat flatter at the origin than in the m — 3 case of 
Fig. |l^, resulting in a slightly less accurate fit to the /3- 
model. However, for the m = 2 case, the obtained simu- 
lated profile is almost indistinguishable from its /3-model 
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radius (h^ ^ Mpc) 

Figure 10. The baryonic number density of tlie cluster (solid 
line) as a function of radius at the time tc = YSSh^^ Myr ago, 
when a photon traversing the centre of the cluster experiences 
the maximum density; m = 3 is assumed. The dashed line is the 
best-fit /3- model (see text). 
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radius (h^ Mpc) 

Figure 11. Log-log plot of the baryonic number density of the 
clusters (solid lines) as a function of radius at the time tc = 
738hg^ Myr ago. Simulations are for m = 2,3 and 4 from right 
to left. The dashed line is the best-fit /3-model for the m = 2 case, 
while the dotted line is for the m = 4 case (see text). 



best-fit from r = up to r ~ IOHq^ Mpc. This can be 
seen in the log-log plot of Fig. |ll], where the density pro- 
files of models m = 2, 3 and 4 are plotted in solid lines 
for O.Ol/io ^ Mpc < r < SOO/Iq ^ Mpc. The dashed line in 
Fig. ^ corresponds to the best-fit /3-model for m = 2. 
In this case, we find no — 0.99 ± 0.15 x /!,,7^''^10'* m~^, 



R 



0.19 ± 0.03/io ^ Mpc and (3 = 0.5 ± 0.05. The best-fit 



/3-model for the m = 4 case is the dotted line of Fig. |ll| 
and corresponds to no — 1.11 ± 0.14 x /Iq^^'^10* m~ , 
R = 0.18 ± 0.03/io ^ Mpc and /3 = 0.67 ± 0.05 which is 
equivalent to the King profile described in equation (^. 

We note that the central density of our simulated clus- 
ter is well defined, unlike some of the {a, /3, 7) profiles of 
equation (Q), which diverge at r = 0. Therefore, in addi- 
tion to the density distribution, we choose to fit the mass 
profile M{r) defined by equation (^ at the time t — tc- 
For both the analytical profiles defined by equation (Q) 
and our simulated cluster, the mass density is defined by 
p{r) = lOhQ^^^nip X n{r), where rUp is the proton mass and 
the factor lOh^^^^ is the ratio of the total mass to the bary- 
onic mass. The mass profile obtained from our model, to- 
gether with the best-fit King and NFW profiles from r = to 
r = 1.5hg^ Mpc, is shown in Fig. ^for the case m = 3. We 
note that our model agrees with masses derived from both 
King and NFW profiles for r > OAhg^ Mpc while a King 
model is favoured for small radii. This is not surprising since 
both the King model and our simulated density profile are fi- 
nite at the origin, whereas the NFW profile diverges at r = 0. 
A least-squares fit of the all the parameters of equation (^ 
gives (q,/3,7) = (1.6,0.67,0.3) which seems to favour the 
simple King profile of equation In this case the num- 
ber density approaches the power law n oc 1/r^ which is in 
agreement with the isoth ermal sphere approximation that 
we made in Section 3.1.2 regarding the /i-dependences. 

It is a particularly satisfying feature of our approach 
that the radial density profile of the perturbation evolves 
into one that matches quite closely those of observed clusters 
and clusters produced in A'^-body simulations. It should be 
emphasised again that the initial conditions assumed a very 
simple finite velocity perturbation with an associated den- 
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radius (h^ ^ Mpc) 

Figure 12. The total gravitational mass of the cluster included 
within a given radius (solid line) at the the time tc = 738/i^^ Myr 
ago, when a photon traversing the centre of the cluster experiences 
the maximum density; m = 3 is assumed. The dashed line is the 
best-fit King profile while the dot-dashed line is the best-fit NFW 
profile. 



sity perturbation. It is very encouraging that these straight- 
forward initial conditions can lead to such realistic density 
profiles as that presented in Figs. 0, and |l|. 

As mentioned above, we define the distance to the clus- 
ter as D = 250/i(^^ Mpc and the cluster core radius to 
be Rc = 0.23^0 ^ Mpc. From our simulations, however, we 
may also measure other important parameters describing 
the cluster at the time when the photon passes through its 
centre. Of particular interest is the turnaround radius Rt, 
which is defined as the radius at which the fluid velocity 
changes from being radially inwards to radially outwards. 
Another important characteristic is the ratio 5 of the clus- 
ter central density to that of the external Universe. We find 
Rt = 34.2/io\ Rt = 16.5/io^ and Rt = 11 Ah^^ Mpc for 
m = 2, 3 and 4 respectively. The density ratio found in all 
three cases is 5 = 7 x 10''. 

We may compare these basic properties with those 
found by Panek (1992) and Quilis et al. (1995). In making 
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Table 1. Previous estimates of cluster properties. D denotes the 
distance to the cluster, Rc is the cluster core radius, Rt is the 
turnaround radius and <5 is the ratio of the central cluster den- 
sity to that of the external Universe. All distances are quoted in 
ftg ^ Mpc 



Author 




D 


Rc 


Rt 


<5 


Panek I 




100 


0.30 


11.7 


674 


Panek II 




100 


1.00 


18.8 


674 


Quilis et al. 




100 


0.23 




3220 


This work (m 


= 2) 


250 


0.23 


34.2 


7000 


This work (m 


= 3) 


250 


0.23 


16.5 


7000 


This work (m 


= 4) 


250 


0.23 


11.4 


7000 



Table 2. The total gravitational mass contained within spheres of 
given radii for a spherical cluster. Mm=2,3,4(f) gives the mass cal- 
culated in this work for m = 2, 3 and 4 models, whereas Mpi(r), 
^PllC'") S'lid A4'Q(r) refer respectively to estimates for the two 
Panek models and the Quilis et al. model. All masses are quoted 
in lO^^hg^ Mq and radii are in h^^ Mpc. 



r 


1 


1.5 


2 


4 




1.45 


2.81 


4.42 


12.6 


Mm=3{r) 


1.01 


1.72 


2.47 


5.84 


Mm=i{r) 


0.78 


1.22 


1.67 


3.55 


Mpi(r) 








1.8 


Mp„{r) 








5.7 


MQ{r) 




0.95 


2.0 


6.2 



such a comparison, however, we must remember that these 
authors specify different initial conditions. The reader is re- 
ferred to their respective papers for details of the assumed 
initial density profiles. In Table ^, we compare the proper- 
ties of the resulting clusters for Panek type I and type II 
models and the Quilis et al. model. 

We see from the table that the turnaround radii for the 
previous model clusters are in reasonable agreement with 
our estimate for fto = 1- However, the density ratio 5 de- 
rived from our simulations in this case is over twice that 
for the Quilis et al. model and over a factor of ten greater 
than those quoted for the two Panek models. This suggests 
that a velocity perturbation of the form assumed here is in 
fact more effective in producing liighly non-linear structures 
than the initial density perturbations assumed by previous 
authors. 

We may also calculate the total mass of the cluster con- 
tained within spheres of various radii and compare the result 
with previous calculations. These are given in Table ^ The 
derived masses for models m = 2, 3 and 4 and agree rea- 
sonably well with Quilis et al. and Panek, particularly for 
m = 3. 

In concluding this section, it must be noted that at later 
times than those illustrated above, the subsequent evolution 
of the cluster is rather unrealistic. As a result of the assump- 
tions of a pressureless fluid and a radial velocity distribution, 
the collapse of the cluster continues unchecked, leading to 
the ultimate formation of a singularity. Hence the formation 
of stable clusters is not admitted by this model. Neverthe- 
less, there are no difficulties in studying photon propagation 



beyond the time the singularity forms, as long as the photon 
has passed the origin by this point. Indeed, it is very useful 
to be able to consider photons that passed through the clus- 
ter shortly before the singularity formed, but which were not 
received by the observer until some time later. The inclusion 
of a baryonic component with pressure will be discussed in 
a forthcoming paper. 



4 CMB ANISOTROPY 

Using the photon propagation equations (Paper I) , we calcu- 
late the predicted CMB anisotropy produced when a photon 
passes through the collapsing cluster. As discussed above, 
however, the assumption of a pressureless fluid and spheri- 
cal symmetry result in the cluster evolution being too rapid. 
Therefore the anisotropies derived in this section should 
properly be considered as upper limits. 

For our present purposes let us ignore primordial 
fluctuations in the CMB and concentrate on secondary 
anisotropies caused by the interaction of CMB photons with 
non-linear structures such as clusters. Three main secondary 
anisotropies can occur; (i) the kinetic Sunyaev-Zel'dovich 
(SZ) effect resulting from the peculiar bulk velocity of a 
cluster relative to the observer, which results in the Doppler 
shift of the CMB blackbody spectrum; (ii) the thermal SZ 
effect produced by the inverse Compton scattering of CMB 
photons to higher energies by hot electrons in the cluster 
gas; (iii) a non-linear gravitational effect on the CMB pho- 
ton as it passes through the cluster. In this section, we will 
be concerned with the last of these, which is often called the 
Rees-Sciama effect (Rees & Sciama 1968). 

A simple, approximate way to understand why such an 
effect should occur is to consider the gravitational potential 
well experienced by the CMB photon. Clearly, if the cluster 
is collapsing as the photon passes through it, then the pho- 
ton must 'climb out' of a well deeper than that into which 
it 'fell', resulting in a net redshift of the photon. However, 
a compensating effect also occurs, since a photon passing 
through the potential well is delayed with respect to a pho- 
ton that did not pass through the cluster. Thus, if these 
two photons arrive simultaneously at the observer, the one 
that passed through the cluster must have been emitted at 
an earlier time, when the Universe was hotter, leading to a 
blueshift of the photon. Thus there is an interplay between 
these two competing effects, and the overall effect can be 
of either sign depending on the details of the evolution of 
the potential well. However, as seen in Paper I (equations 33 
and 36), the gauge theory approach offers a clear alternative 
understanding of the physical processes involved in the Rees- 
Sciama effect, showing that the microwave decrement is di- 
rectly related to the velocity gradient du{t,r)/dr encoun- 
tered by the photon rather than the potential well variation. 
This raises questions on issues related to the generalisation 
of available models to the non-spherical case and models in- 
cluding a pressure component as discussed in Quilis & Saez 
(1998). 

Fig. ^ shows the gravitational CMB anisotropy due to 
the cluster described in the last section, for m = 3. The 
anisotropy is plotted in AT/T as a function of projected 
angle 6 on the sky from the centre of the cluster. The max- 
imum anisotropy occurs at the centre of the cluster and 
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Table 3. Various estimates of the anisotropy AT/T compared 
with our three estimations for models m = 2, 3 and 4. 



6 (degrees) 

Figure 13. The predicted CMB anisotropy (solid line) due to 
the gravitational effect of the cluster described in Section ^, for 
m = 3. The anisotropy is plotted in units of AT/T as a function 
of the projected angle on the sky from the centre of the cluster. 
The dashed line shows the fluid velocity in arbitrary units at the 
point of maximum density along any given photon trajectory. 



-800 -600 -400 

time (h(j ^ 10^ years) 



-1000 -BOO -600 -400 -200 

time ^ 10^ years) 

Figure 14. Top: the change in AT/T for a CMB photon passing 
through the centre of the cluster, as seen by an observer comoving 
with the fluid at any given time for the m = 3 model. Bottom: 
the same curve plotted between smaller limits. AT for observers 
inside the cluster is defined with respect to the temperature of the 
external Universe at the same time. The time t = corresponds 
to the current epoch. 



Authors 




AT/T 


Panek type I 




-0.15 X lO-'^ 


Panek type II 




-0.6 X 10~5 


Quilis et al. 




-1.2 X 10^5 


Chodorowski I 




-0.26 X 10""' 


Chodorowski II 




-0.77 X 10-"' 


Nottale 




« -10 X 10""' 


This work (m = 


2) 


-5.8 X 10^5 


This work (m = 


3) 


-0.93 X 10-5 


This work (m = 


4) 


-0.4 X 10-5 



has the value AT/T = -0.93 x 10"^ -phe decrement de- 
creases rapidly with projected angle on the sky and reaches 
zero at 9 = 4.4 degrees. For larger angles the anisotropy 
becomes slightly positive indicating a net blueshift of the 
CMB photons. This positive effect reaches a maximum of 
AT/T = 0.25 X 10"^ at e» = 6.8 degrees before slowly tend- 
ing towards zero as 9 increases. We note that, using the 
/i-dependences stated in Section 3.1.2, the obtained AT/T 
is independent of ho. Also plotted in Fig. |l3|is the velocity 
in arbitrary units at the point of maximum density along 
any given photon trajectory. It is interesting to note that 
the crossover between the negative CMB anisotropy and the 
positive one at 6' = 5 degrees occurs close to the point at 
which the velocity changes from being radially inwards to 
radially outwards. 

Anisotropics of a similar shape are also found for the 
m — 2 and 4 cases and we find the maximum decrements 
to be AT/T = -5.8 x 10"^ and AT/T = -0.4 x 10"^ 
respectively. The result for m = 2 is somewhat larger than 
that found for m — 3, which is to be expected since the 
cluster obtained for m = 2 is larger even though it sa tisfies 
the same observational requirements stated in Section 3.1.2. 



We can compare the central decrements calculated 
above with those of previous authors. Results are reported 
in Table ^ for Panek (1992) type I and type II clusters, 
Quilis et al. (1995), Chodorowski (1991) models I and II 
and finally for Nottale (1984), who used the Swiss Cheese 
model to predict a considerably larger central decrement of 
AT/T ~ — 10 X 10-5. However, this last value corresponds 
to very dense, rapidly collapsing objects and the resulting 
cluster mass is much greater than those observed. 

We see that our predicted central decrement for m — 3 
is in rough agreement with that quoted by Quilis et al., 
which is slightly larger than those found by Panek and 
Chodorowski. However, as discussed above, the effect pre- 
dicted here has to be regarded as an upper limit, particularly 
in our m = 2 model where AT/T = -5.8 x lO'^. 

In Fig. |l^, we plot the change in AT/T of a CMB pho- 
ton passing through the centre of the cluster, as seen by an 
observer comoving with the fluid at any given time, for the 
m — 3 model. The temperature difference is measured rela- 
tive to CMB photons that have not interacted with a clus- 
ter. The top plot illustrates the wide range of anisotropy 
produced by the interaction with the cluster. The bottom 
plot displays the same curve plotted between smaller limits 
and shows the resulting AT/T = —0.93 x 10"^ as discussed 
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above. We note that AT/T ~ at the origin which corre- 
sponds to tc — 738/i^^ Myr ago. This plot also illustrates 
that AT/T does not change when the photon is in the region 
exterior to the cluster. Therefore, as long as the observer lies 
in the external Universe (i.e. sufficiently far away from the 
cluster centre), the magnitude of the Rees-Sciama effect is 
independent of the distance to the cluster. In the case where 
the observer lies sufficiently close to the cluster that he is 
decoupled from the Hubble flow, predictions of the Rees- 
Sciama effect magnitude could be highly inaccurate, as seen 
from the large gradients of Fig. |l^ (bottom panel). In the 
case where m = 2, the sharp gradient in AT/T that occurs 
at t ~ 500/i^^ Myr ago in Fig. |l4| is shifted to a time close 
to the present day. Indeed, this is the reason why we had to 
place our clu ster a t an effective distance of _D = 250/1,^^ Mpc 
(see Section 3.1.2). 



5 APPLICATION TO HIGH-REDSHIFT 
CLUSTERS 

We have so far concentrated on the formation of rich, low- 
redshift Abell clusters in order to compare our results with 
previous authors. We have commented, however, that such 
clusters are most likely virialised to some extent and there- 
fore are not well modelled by a spherical free-fall collapse. 
Nevertheless, such a model may provide a better descrip- 
tion of the dynamical state of high-redshift clusters. We 
would expect clusters at z ^ 1 still to be in the process of 
formation, and therefore may be reasonably approximated 
as quasi-spherical objects with large radial infall velocities. 
Massive clusters at high redshift have recently been identi- 
fied, either by X-ray measurements, gravitational lensing or 
spectroscopic identification {z = 0.83, Luppino and Kaiser 
1997; z = 0.996, Deltorn et al. 1997). 

The existence and distribution of high redshift clus- 
ters are of great importance regarding the discrimination 
between different cosmological scenarios. The observational 
aspect of the problem has developed only recently. Indeed 
identifying bound structures at high redshift is a real ob- 
servational challenge not only because of the decrease in 
brightness with distance but also because of the high degree 
of confusion due to the projection of numerous foreground 
and background objects. As a search strategy, Jones et al. 
(1997) proposed to use the Sunyaev-Zel'dovich (SZ) effect 
to detect the presence of high-redshift clusters, the main 
argument being that the magnitude of the SZ effect is in- 
dependent of redshift. Since the Rees-Sciama and the SZ 
effects usually add up to form a total observed CMB decre- 
ment, it is essential to know precisely the magnitude of the 
Rees-Sciama effect so that the SZ contribution can be eval- 
uated accurately. Therefore, a model such as ours, which 
allows us to understand the details of the Rees-Sciama ef- 
fect, should prove useful in the context of programmes such 
as that carried out by Jones et al. (1997) or similar future 
work. 

In order to apply our model to some on-going observa- 
tions, we choose to investigate the properties of high-redshift 
clusters which may model the microwave decrement reported 
towards the quasar pair PC1643-I-4631 A&B (Jones et al. 
1997). This decrement is observed at 15 GHz and lies be- 
tween a pair of quasars at redshifts z = 3.79 and z — 3.83 















S 











4" 



Figure 15. Geometrical arrangement of a source-lens system. 
S denotes the plane of the source, which is comoving with the 
Hubble flow. A photon is emitted from an arbitrary position on 
this plane towards the observer O. The emission point is defined 
by the observed angle (in presence of the lens) or by the observed 
angle /3 (in absence of the lens) and a = 9 — (3. Dg^ denotes the 
angular-size distance from the point x to the point y. 



separated by 198 arcsec on the sky. Jones et al. assume the 
decrement to be a thermal Sunyaev-Zel'dovich effect due 
to an intervening cluster and estimate the minimum cen- 
tral temperature decrement to be AT/T = —2.1 x 10"*. 
However, no cluster is evident in either ROSAT X-ray ob- 
servations or infrared R-, J- and K-band images taken by 
the WHT and the UKIRT (Saunders et al. 1997). The mag- 
nitude limits obtained lead the authors to suggest that any 
cluster causing the decrement must lie at z ^ 1. 

Modelling of the SZ observations, assuming ho = 0.5, 
suggests that the putative cluster has a core radius of 
Rc ~ 0.3 — 0.5 Mpc and a gas mass inside a 2- Mpc radius of 
2.0-3.5x10^* Mq. This implies a total mass, including both 
luminous and dark matter, greater than 10^* Mq. Saunders 
et al. further point out that for such a massive cluster at 
= 1 — 2, the Einstein ring radius for sources at z « 3.8 is 
approximately 100 arcsec. Indeed, simple gravitational lens 
modelling can simultaneously produce the observed SZ ef- 
fect and make the true positions of the two quasars virtually 
coincident. The suggestion that the observed quasar pair are 
in fact images of the same object is supported by the remark- 
able similarity between the quasar spectra, apart from their 
one-percent redshift difference. Saunders et al. suggest that 
this difference could be explained in terms of quasar evolu- 
tion over a delay between the two lightpaths of ~ 10^ yrs. A 
further similar observation has been reported by Richards 
et al. (1997) where a microwave decrement is also detected 
towards a pair of quasars. 

Using the model described in Section ^, we now con- 
sider the formation of a massive cluster located ai z — 1 
and investigate both the gravitational lensing effects and 
the microwave decrement produced by such a cluster 



5.1 Gravitational lensing 

We calculate the evolution of the mass distribution and the 
photon paths as described in Sections ^ and ^ using the fol- 
lowing parameters: Q{t) = 1; H{to) — 50 km s"'^ Mpc~^; 
the maximum baryonic number density along the photon 
path is taken to be 10* m~^; the total gravitational mass 
over baryonic mass ratio is fixed to 10; the initial perturba- 
tion is for m = 3; the observed cluster redshift is taken to be 
z = 1 and the core radius is 0.45 Mpc (52 arcsec). The source 
of photons corresponding to the observed pair of quasars is 
comoving with the Hubble flow and lies at a redshift z = 3.8. 
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8 (arcsec) 

Figure 16. The reduced deflection angle a as a function of view- 
ing angle 9. The cluster lies at z = 1 and the source at z = 3.8. 



10 100 1000 10* 

9 (arcsec) 

Figure 17. The inverse of the magnification factor as a function 
of the viewing angle 8. The magnification becomes infinite at 
fie ~ 46" and at ~ 100". 



We note that the assumptions made here concerning the size 
and densit y of the cluster are somewhat different than in Sec- 
tion IS.l.l In order to fit PC1643-I-4631 A&B observations, 



we find it necessary to model such a rich cluster. 

We now consider the simple case where a (slightly ex- 
tended) source lies with its centre on the same line of sight as 
the cluster. For a spherically symmetric cluster, the resulting 
image is radially symmetric and any position in the image 
is defined by the radial angle 6 measured from the cluster 
centre. We can define the reduced deflection angle a{6) by 
using the standard lens equation (Refsdal, 1964; Blandford 
& Narayan, 1992): 



aie) 



-m, 



(9) 



where 9 and /3 are the observed viewing angles of the same 
point on the source plane with and without the lensing mass 
respectively. We note that in the context of our model de- 
scribed in Section |^, 9 and /3 are the relevant values of the 
angle Xi as observed by the comoving observer O of Fig. |l5[ 
In the geometrical arrangement of Fig. the reduced deflec- 
tion angle can be read as the convergence (i.e. radial stretch- 
ing) of the source. For various image positions 9, we compute 
the corresponding angle /3 and so the convergence a{9) in 
the case of our evolving gravitational lens. We find that the 
lensing effect occurs up to rather large angles ( ~ 3 degrees) , 
as seen in Fig. |l6| . 

Another interesting characteristic of our source-lens sys- 
tem is its magnification factor fi{9) defined to be the ratio 
between the observed solid angles of the same unit area at 
the source plane, with and without the lens. Using (^ in the 
circularly symmetric lens case, we find 

N_da(S)n 

(10) 

Using Fig. ^ we may compute the corresponding magiufi- 
cation, which is displayed in Fig. |l^. We can see from ( p^ 
that there are two cases where the magnification diverges. 
First, on the Einstein ring radius 6^, we have 9e — a(Se) = 
/3{9e) = 0. Any photon observed with 9 — 9^ originates from 
the same emission point at the centre of the source plane. 
Therefore the solid angle observed in the absence of the lens 
shrinks to zero and the magnification becomes infinite. Sec- 
ond, the magnification also diverges when da/d6 — 1 (in 




(arcsec) 

Figure 18. Lensing image of four sources A, B, C & D located at 
the same redshift z = 3.8. Aq, Bq, Cq, Dq denotes the positions 
of the sources as observed in the absence of the cluster. Ai, A2, 
Bi, Ci and Di are the corresponding lensed images. 



this case we have d/3/d6 = 0); we denote the angle at which 
this occurs by 9c. 

On the source plane, the caustic defined by the first 
case is a critical point whereas the caustic defined by the 
second is a ring. These caustics are called the point caustic 
and outer caustic, respectively. A point source at the point 
caustic produces a ring of radius 9e on the image plane (i.e. 
the Einstein ring); a point source on the outer caustic is 
observed as a double source with one radially elongated and 
reversed image &t9 = 9c and a second tangentially elongated 
image at > 6'e. In general a point source located outside 
the outer caustic produces a unique lensed image while a 
triple image is obtained if the source lies within the outer 
caustic (Blandford & Narayan, 1992). 

As an illustration of the lensing properties of our clus- 
ter, in Fig. ^ we show a lensed image of four extended 
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0.6 O.B 1 1.2 1.4 

Pmax (10* P m"') 

Figure 19. Einstein ring radius contours over a range of pmax 
and Rc values for, from top to bottom, 0e = 150", 9e = 100" 
and 6^ = 50". The solid lines denotes the contours computed 
assuming a simple static lens with an King density profile. The 
points are exact calculations for an evolving lens (squares for 6e = 
150", circles for Be = 100" and triangles for Be = 50"). 



Table 4. Properties of clusters in Fig. [l9[ Be denotes the Einstein 
ring radius; Te is the electron temperature in the cluster; ATgz is 
the Sunyaev-Zel'dovich decrement and ATag is the Rces-Sciama 
decrement. 



Cluster 


9e 


Te (K) 


ATsz (^iK) 


ATrs (pK) 


A 


150" 


2.12 X lO'' 


523 


306 


B 


150" 


1.54 X 10^ 


466 


283 


C 


150" 


1.36 X lO'^ 


459 


274 


D 


100" 


2.27 X 10^ 


491 


201 


E 


100" 


1.95 X 10^ 


464 


185 


F 


100" 


1.78 X lO'' 


458 


178 


G 


100" 


1.63 X 10^ 


459 


172 


H 


100" 


1.56 X 10^ 


463 


169 


I 


50" 


2.64 X lO'' 


464 


117 


J 


50" 


2.26 X 10^ 


466 


87 


K 


50" 


2.12 X lO'' 


490 


79 



sources A, B, C & D. Each source is located at the same 
redshift z = 3.8 and is assumed to be circularly symmetric. 
Sources B, C, D lie outside the outer caustic whereas A is 
positioned on the caustic itself. We note the double image 
(almost triple on the figure) produced by A and its tangen- 
tial/radial characteristics. The positions of the sources as 
observed without the lens are indicated so that the lensing 
effect of the cluster can be fully appreciated. 

It is also interesting to note the effect of the lens on a 
background of sources at a given redshift with a projected 
number density (in the absence of the lens) of No{8). In the 
presence of the lens the observed number density becomes 
N(e) = No{e)/^i, and from Fig. |l7|we see that N(e) < No{e) 
for 0^1 degree. 



5.2 Microwave decrement 

In this section we investigate the contribution of the Rees- 
Sciama effect to the total temperature decrement reported 



towards the quasar pair PC1643+4631 A&B and discuss the 
possibility that the quasars are gravitationally lensed images 
of the same object (Dabrowski et al. 1997; Dabrowski 1997). 

In order to obtain any quantitative results we need 
to constrain our model parameters. We suppose that both 
quasars have a redshift z — 3.8 and that the cluster lies 

we assume that Q{t) 



at z 
H{to 



1. As for Section 5.1 



1, 



50 km s~^ Mpc"'^ and that the total gravitational 
mass over baryonic mass ratio is 10. The initial perturbation 
is for m = 3. Two parameters remains to be fixed: pmax, the 
maximum baryonic density encountered by an observed pho- 
ton and the characteristic core radius Rc of the cluster at the 
time the photon reaches the point of maximum density. To 
fix these parameters, we have two observational constraints: 
the quasar pair separation of 198 arcsec and the total flux 
observed by the Ryle Telescope of —380 /xJy (Jones et al. 
1997). 

We first consider the lensing separation. Naturally, our 
spherically-symmetric model cannot account for multiple 
images that do not lie on the same radial line from the cen- 
tre of the cluster, as in the case of PC1643-I-4631 A&B. We 
note here that a more suitable elliptical setup might present 
a longer projected line of sight and so enhance the CMB 
decrement effect. For the purpose of this paper, we consider 
the Einstein diameter 2x9sto be characteristic of the lensing 
power of the cluster and compare this with the separation 
of the quasar images. For each pair of values pmax and Rc, 
which describe the cluster central density and core radius, we 
obtain a value for the Einstein ring radius 9e. The solid lines 
in Fig. |l^ are contours in (pmax, -Rc)-space corresponding to 
Einstein ring radii of 6c — 50, 100 and 150 arcsec respec- 
tively, calculated without using our model but assuming a 
static lens with a density distribution described by a King 
profile (see equation |^ . In this case the Einstein ring radius 
may be found simply by solving numerically the standard 
equation 



(11) 



where M{9c) is the projected mass included within the ra- 
dius 6c, Dg'^ denotes the angular-size distance from the ob- 
server to the cluster, Dg" the angular-size distance from the 
observer to the source and Dg" the angular-size distance 
from the cluster to the source (see Fig. |l5| ). These distances 
were computed using the formulae given by Blandford & 
Narayan (1992). 

In Fig. I19L we also plot 11 selected cluster setups 
that have been computed using exact calculations from our 
evolving spherical cluster model. Since these points repre- 
sent evolving lenses, and the corresponding clusters do not 
have exact King density profiles, we see that these points 
lie slightly away from the contours. The points plotted 
have been chosen so that A, B, C denote lenses for which 
Oc = 150"; D, E, F, G, H have Oc = 100" and I, J, K have 
Oc = 50". 

Secondly, we consider the total fiux observed by the 
Ryle telescope between the pair of quasars. Table ^ shows 
the electron temperatures Tc required in order that the clus- 
ters in Fig. |l^ each produce a total observed SZ flux of 
—380 /iJy. Also listed in the table are the corresponding SZ 
and Rees-Sciama temperature decrements for each cluster. 
We note that the ratio of the Rees-Sciama effect as com- 
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Figure 20. Left: a realisation of CMB fluctuations in a standard inflationary CDM model. Right: the same realisation but with a rich 
cluster at its centre, producing a Rees-Sciama temperature decrement and gravitational lensing effects. 



pared to the SZ effect varies from 0.61 for cluster B to 0.16 
for cluster K. Thus in most cases the Rees-Sciama effect is 
not negUgible and might contribute significantly to the total 
observed CMB decrement. As pointed out in the Introduc- 
tion, however, we note that these values are to be taken as 
upper limits. 

We take cluster F as a typical example of a cluster which 
may explain the observations of Jones et al. (1997). The total 
mass contained within spheres of various radii for cluster F 
are: 7.8 x 10^^ Mq within 2 Mpc; 1.9 x lO^** Mq within 
4 Mpc; and 4.6 x 10^*^ M© within 8 Mpc. The ratio of the 
central density to the density of the external Universe at the 
time the photon passed through the cluster is 5 = 4440 and 
the turnaround radius Rt — 25.3 Mpc. 

We may also consider the possible time delay between 
the light paths for the two quasar images PC1643-I-4631 
A&B. A reasonable approximation may be obtained in our 
model by considering a source located directly behind the 
cluster at z = 3.8 and calculating the time delay between 
a photon from the source travelling straight through the 
centre of the cluster and one which follows a lensed path, 
appearing at the Einstein ring radius 6e- We note that our 
model takes full account of all the relativistic effects caused 
by the motions of the source and observer in the Hubble 
flow and the gravitational effects of the collapsing cluster. 
We find that the proper time delay in the frame of the source 
is ~ 150 years. This is rather a brief period for the quasar 
to evolve sufficiently to explain the small redshift difference 
of the two quasar images. 

Finally, we note that observations of the PC1643-I-4631 
field are still in progress. For example, Haynes et al. (1998) 
carried out deep optical imaging with the William Herschel 
Telescope and show that the excess of faint blue galax- 
ies in the field might be consistent with a cluster lying at 
2: ~ 2. Furthermore, very recent ROSAT observations may 
suggest that any intervening cluster producing the tempera- 



ture decrement towards the quasar pair PC1643-1-4631 A&B 
should in fact have a redshift z > 3 (Kneissl 1997; Kneissl, 
Sunyaev and White 1998). We have therefore repeated our 
analysis for a cluster lying at a redshift greater than 1 and 
chosen the (maybe conservative) value of z = 2. In this case 
we consider a cluster with a baryonic number density of 
10* m~^ and a core radius Rc = 0.74 Mpc; this ensures an 
Einstein ring radius of 100 arcsec. The total mass contained 
within a sphere of 2 Mpc is 1.74 x 10^® Mq, the density ratio 
5 — 1313 and the turnaround radius Rt — 22 Mpc. In order 
to retrieve the observed SZ flux of —380 /iJy, we require the 
electron gas temperature of the cluster to be 1.3 x 10^ K. 
This gives an SZ temperature decrement of 547 fiK, while 
the upper limit of the Rees-Sciama temperature decrement 
is in this case ~ 600 fiK (i.e. it is of the same order as the 
thermal SZ effect). For such a cluster, we find that the lens- 
ing time delay, as described in the previous paragraph, is 
approximately 230 years in the frame of the emitting source 
placed at a redshift z = 3.8. 



6 GRAVITATIONAL LENSING OF 

PRIMORDIAL CMB ANISOTROPIES 

We discussed in the previous section a possible unified model 
for the formation of a high-redshift cluster with the proper- 
ties required to produce a microwave decrement and gravi- 
tational lensing effects consistent with the observations to- 
wards the quasar pair PC1643-I-4631 A&B. It is also of some 
interest, however, to investigate the effects of such a cluster 
on primordial CMB fluctuations. 

In this Section, we therefore consider the non-linear 
gravitational effects on a typical spectrum of primordial fluc- 
tuations of the massive cluster F, discussed in the last sec- 
tion, at a redshift z — 1. These effects include both the 
gravitational lensing and shift in energy of CMB photons 
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passing through the cluster. We do not, however, include 
any SZ effect that the cluster may produce. 

The effect on primordial fluctuations is investigated by 
propagating through the evolving cluster a population of 
CMB photons due to a typical primordial CMB field. It 
will be seen from Fig. ^ that the reduced defiection angle 
produced by this cluster falls to zero only at large angles of 
about 3 degrees. We therefore take as our typical primordial 
CMB field a 4 X 4-degree realisation of CMB fiuctuations 
in a standard infiationary Universe dominated by cold dark 
matter. 

Fig. ^ shows that effect of the cluster lying at the cen- 
tre of the field of primordial CMB fluctuations. We see that 
the cluster produces a pronounced Rees-Sciama decrement 
as well as causing a slight radial stretching of the CMB fluc- 
tuations which extends as far as the edge of the field. 

6.1 Massive cluster abundances 

We may also consider the effect of a population of such clus- 
ters on the CMB power spectrum. As seen in previous work 
(e.g. Eke, Cole & Frenk 1996) standard models using the 
Press-Schechter formalism in a S7o = 1 Universe predict 
no structure of mass ~ 10^^ at redshift greater than 
z ~ 0.8. However counts of clusters are sensitive to Qq so 
that open cosmological models with Qo ~ 0.3 allow the pres- 
ence of massive clusters at 2: > 1. Eke et al. (1996) extended 
the Press-Schechter formalism to flat cosmological models 
with a cosmological constant A so that Qa + fio = 1, where 
Qa = A/{3Hi). For a given value of f2o, the redshift distri- 
bution of massive clusters is similar in both open and flat 
cosmological models. We assume here Qa + Qo ~ ^ and a 
cold dark matter model. We have used the formalism and 
normalisation described in Eke et al. (1996). The number 
density of clusters with mass between M and M + dM at a 
given redshift z is given by n{z, M)dM with 



n{z, M) 



\MsJ 



poaSc{z ) 

(78 M2 



■ exp 



L "8 V J«s ; 



(12) 



where Mg is the mass contained in spheres of radius 



8/1^^ Mpc; erg = 0.52r2, 



-0.52+0.1300 



is the rms linear fluc- 



tuation amplitude within 8/ig ^ Mpc spheres; 
3ii-| 



po 



SttG 



(I-^^a) 



(13) 



+ 



is the present mean density of the Universe; a = (0 
0.4r)/3 (see White, Efstathiou and Frenk 1993) for F = 
f2o/io = 0.25; and 5c{z) is the density threshold for collapse 
as derived in Eke et al. (1996) in the case Q.a + Qo = 1. 
The number of clusters of mass from M to M -I- dM within 
a spherical shell extending from z to z + dz is given by 
n{z,M)dMdV with dV = AirD^dr where D is the effective 
distance and dr is the comoving radial coordinate element 
equal to —(1 + z)cdt. Therefore 

dt 



dV 



-^AncD^ (1 + z] 



dz 



dz 



(14) 



where dt/dz and D are derived from Fukugita et al. (1992) 
for the Q.A -I- fio = 1 case: 



dt_ _ 

dz 

and 



Ho{l + z)^?,noz(l + z) + l + ftoz'-i 



(15) 




Figure 21. The number of clusters with mass larger than 
10^^ h^^ Mq per unit redshift in the whole sky, as predicted by 
the Press-Schechter formalism. We assume Qq = 0.3, Slo-f Ha = 1 
and the normalisation of Eke et al. (1996) 




Figure 22. The number of clusters with mass larger than 
10^^ h'^^ Mq at a distance further than the given redshift z in 
the whole sky, as predicted by the Press-Schechter formalism. 
The lower, middle and upper curves are for Qg = 0.3, Qg = 0.23 
and Qq = 0.17 respectively and fig + Qa = 1 is assumed in each 
cases. The normalisation is that of Eke et al. (1996) 



dz 



^ Ho.jg + zf' + 



(16) 



We use the above formulae and find the distributions 
given in Figs. |2l] & ||. 

Fig. pil shows that for fio ~ 0.3 only very few clusters 
more massive than 10^^ Mq are expected at redshift larger 
than 1. Fig. ^ gives the estimated number of lO^^M© clus- 
ters in the whole sky as a function of redshift, for various 
value of fio- For Q.q — 0.3, 0.23 and 0.17 there are, in the 
whole sky, ~ 4, 40 and 400 clusters respectively at 2 > 1 and 
more massive than IO^^A^q. We note that the presence of 
massive objects as distant as PC1643 may suggest that re- 
sults using Press-Schechter formalism in a standard cosmol- 
ogy underestimate massive cluster counts at high redshift. 
This is the reason why, in addition to the realistic assump- 
tion Qo = 0.3 (e.g. Webster et al. 1998), we are considering 
here Qo ~ 0.23 and ilo ~ 0.17 which give abundances respec- 
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tively 10 times and 100 times higher than those expected for 
^0 = 0.3. 



6.2 Effect on the CMB power spectrum 

The presence of 4 clusters in the whole sky corresponds 
roughly to one cluster observed in every 100 x 100-degree 
field on the sky, while 40 clusters correspond to one clus- 
ter every 31 x 31-degree field and 400 to one cluster ev- 
ery 10 X 10-degree field. We performed, for these 3 cases, 
Monte-Carlo simulations of the gravitational effect of the 
clusters on realisations of primordial CMB fluctuations, and 
measured the corresponding power spectra. Averaging these 
measured spectra, we obtain an estimate of the ensemble av- 
erage power spectrum of CMB fluctuations in the presence 
of each proposed cluster population. The results are shown 
in Fig. p3l where they are compared to the ensemble aver- 
age power spectrum predicted in the absence of clusters. We 
note that our cluster formation model has been developed 
for flo = 1 and Qa ~ 0, while abundances have been esti- 
mated for fio + = 1 (see Section 6.1). We will investigate 



the effect of on our model of Q,o ^ 1 and Q,a 7^ in future 
work. 

It is clear that the effect is negligible if only 4 clus- 
ters are present in the sky. However, for abundances where 
40 or 400 clusters are present in the whole sky, which corre- 
spond to r^o ~ 0.2, the effect is more pronounced and results 
in smoothing out of the Doppler peaks in the inflationary 
power spectrum. Therefore, the possibility of such an effect 
should be taken into account when determining cosmologi- 
cal parameters from future CMB observations. These results 
are in reasonable agreement with the weak lensing analytical 
calculations of Seljak (1996) or Martmez-Gonzalez & Sanz 
(1997). 



7 CONCLUSIONS 

We apply a new model for the formation of nonlinear cos- 
mic structures (Paper I) to the collapse of spherical clusters 
of galaxies. The external, expanding Universe and the col- 
lapsing cluster are governed by the same pressureless fluid 
equations. These equations are exact general-relativistic so- 
lutions of Einstein equations and no approximations have 
been made. The initial conditions for the fluid at an early 
epoch (z = 1000) are very simple: we impose a finite pertur- 
bation on the fluid velocity field that is determined by only 
three parameters and the corresponding density perturba- 
tion is inferred assuming that the perturbation arose from 
primordial fiuctuations. 

In Section |^ we studied the formation of a cluster at a 
redshift z = 0.08 and found that density profile and mass 
distribution of the resulting cluster are realistic. We also 
computed several characteristic quantities for the cluster, 
such as the total mass contained within spheres of various 
radii; the ratio of the central density to that of the external 
Universe {5); and the turnaround radius (Rt)- Comparing 
our results with previous authors (e.g. Panek 1992; Quillis 
1995), we find reasonable agreement. 

Since photon paths are also easily calculated in our 
model, in Section ^ we studied the gravitational effect of the 




Figure 23. Effect on the CMB power spectrum of a population 
of rich cluster of galaxies. Upper panel; the solid curve is the 
unperturbed power spectrum while the dashed curve represents 
the 'lensed' spectrum if 400 clusters are present in the whole sky. 
Lower panel: the curves are the difference between the predicted 
unperturbed power spectrum and the 'lensed' spectrum for three 
different cluster abundances corresponding to 400 (dashed curve), 
40 (dashed-dot curve) and 4 (dotted curve) clusters in the whole 
sky. 



collapsing cluster on CMB photons (i.e. the Rees-Sciama ef- 
fect). For a photon traversing the centre of the cluster, we 
found a central temperature decrement AT/T ~ — 1 x 10~^ 
which is in reasonable agreement with previous estimates. 

Since our model is most relevant to clusters with large 
infall velocities, in Section ^ we apply it to clusters with 
a redshift oi z > 1. Indeed, such high-redshift clusters are 
more likely to be in a state of free-fall collapse than the low- 
redshift clusters considered in Sections ^ and ^. In particular, 
we use our model in an attempt to describe the microwave 
decrement reported towards the quasar pair PC1643-I-4631 
A&B (Jones et al. 1997). Since the quasar pair is possi bly 



lensed (Saunders et al. 1997), we investigated in Section 5.1 
the lensing properties of a cluster which may explain the ob- 
servations. We find that for such a cluster lensing occurs out 
to large projected angles 9 from its centre, with an appre- 
ciable effect still visible at S = 2 or 3 degrees (see Fig, 



In Section 5.2 we consider the relative contributions of 



the Rees-Sciama and thermal SZ effects to the microwave 
decrement observed towards PC1643-I-4631 A&B, and show 
that the Rees-Sciama effect might contribute significantly 
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for clusters that can simultaneously produce the required 
lensing properties discussed by Saunders et al. (1997). At 
z = 1, such a cluster would have a typical central number 
density of 10* m~^ and a core radius of Rc = 0.45 Mpc 
(52 arcsec). The total mass contained within a sphere of ra- 
dius 2 Mpc is then 7.8 x lO^'^ M©. We also find, however, that 
in this scenario the time delay between the light paths for the 
two quasars images PC1643-I-4631 A&B is approximately 
150 years, as measured in the frame of the (lensed) quasar. 
This period might be rather short to explain the slight red- 
shift difference between the two quasar images. Following 
very recent ROSAT observations of PC1643-I-4631 A&B by 
Kneissl (1997), which suggest that any intervening cluster 
should be at an even greater redshift, we also repeat our 
calculations for a similar cluster at a redshift z — 2. We find 
a typical core radius of Rc = 0.74 Mpc (90 arcsec) for a cen- 
tral number density of 10* m~'^. The total mass contained 
within a sphere of radius 2 Mpc is 1.74 x 10^*^ Mq. 

Finally, in Section ^, we consider the effect on primor- 
dial microwave background fluctuations of a population of 
massive clusters, such as that described in Section ^ We 
find that in the case of cluster abundances corresponding to 
a fio ~ 0.2 and fio+^^A = 1 cosmological model, the Doppler 
peaks of the CMB power spectrum are slightly smoothed out 
by the lensing effects (see Fig. ^3|), confirming weak lensing 
the results in Seljak (1996) and suggesting that this effect 
should be taken into account when determining cosmological 
parameters. 
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